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1 Introduction 



The discovery of an approximate SU{3) flavour symmetry of strong interactions PlEl predates 
tlie quark model and has been of paramount importance in all subsequent developments. This 
symmetry has its origin in the smallness of up, down and strange quark masses with respect 
to the QCD scale Aqcd- In static hadron properties, such as masses, magnetic moments, 
decay constants, etc., it is accurate to about 1 to 3% for quantities related by isospin and 
to about 20% for those related by U- and V-spin. The breaking of SU{3) in dynamical 
observables can be larger and up to now is not fully understood. One particularly striking 
example is the weak radiative decay S — >■ J97: the experimental azimuthal asymmetry in this 
decay is = —0.76 ± 0.08 j3j, which according to Hara's theorem |lj implies 100% SU{3) 
symmetry violation. 

In recent years S'?7(3)-symmetry-breaking effects in heavy-meson decays have attracted 
increasing interest. These processes can be treated in heavy-quark expansion, which has 
proved a very powerful theoretical tool, so that in some cases, for instance weak radiative 
decays, B ^ p'j vs. B —>■ K*^, the uncertainty in SU{3) breaking is actually the dominant 
source of theoretical error. The particular challenge of such processes is that, in the presence 
of a hard scale, hard exclusive reactions are dominated by rare configurations of the hadrons' 
constituents: either only valence-quark configurations contribute and all quarks have small 
transverse separation (hard mechanism), or one of the partons carries most of the hadron 
momentum (soft or Feynman mechanism) The size of S'f/(3)-breaking effects in such rare 
configurations cannot be deduced from the symmetry breaking in static properties, where 
the bulk of the wave-function contributes. 

Hard contributions are simpler to treat than their soft counterparts and their structure 
is well understood, see Ref. jEj. They can be calculated in terms of the hadron distribution 
amplitudes (DAs) which describe the momentum-fraction distribution of partons at zero 
transverse separation in a particular Fock state, with a fixed number of constituents. DAs 
are ordered by increasing twist; the leading-twist-2 meson DA 02;p, which describes the 
momentum distribution of the valence quarks in the meson P, is related to the meson's 
Bethe-Salpeter wave function (pp^BS by an integral over transverse momenta: 



Here u is the quark momentum fraction, Z2 is the renormalization factor (in the light-cone 
gauge) for the quark-field operators in the wave function, and fi denotes the renormalization 
scale. The study of the leading-twist DA of the pion has attracted much attention in the 
literature, whereas the status of SU (3)-breaking effects that are responsible for the difference 
between the kaon and the pion DAs has been controversial for a while . These corrections 
have been recently reconsidered in the framework of QCD sum rules in Refs. jHl El ^1 ITT| . 
with a consistent picture finally emerging. We will give a short review of these developments 
below. 

Higher-twist DAs are much more numerous and describe either contributions of "bad" 
components in the wave function, or contributions of transverse motion of quarks (antiquarks) 
in the leading-twist components, or contributions of higher Fock states with additional glu- 
ons and/or quark-antiquark pairs. Within the hard-rescattering picture, the corresponding 
contributions to the hard exclusive reactions are suppressed by a power (or powers) of the 
large momentum Q, which explains why they have received less attention. 
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In turn, soft contributions are intrinsically non-perturbative and cannot be further re- 
duced (or factorized) in terms of simpler quantities without additional assumptions. At 
present, they can only be estimated using light-cone sum rules [IS UHl Clj , see Refs. [T5| IT6] 
for sample applications to heavy quark decays. In this technique soft contributions are ex- 
tracted from the dispersion relations for suitable correlation functions, by introducing an 
auxiliary "semi-hard" scale (Borel parameter) at which the two different representations of 
the correlation function, in terms of quarks and in terms of hadronic states, are matched. 
In calculations of this kind, the necessary non-perturbative input again reduces to DAs, and 
the higher-twist DAs play a very important role, since they are only suppressed by powers 
of the Borel parameter and not by powers of the hard scale Q. The crucial point and main 
technical difficulty in the construction of higher-twist DAs is the necessity to satisfy the exact 
equations of motion (EOM), which yield relations between physical effects of different ori- 
gins: for example, using EOM, the contributions of orbital angular momentum in the valence 
component of the wave function can be expressed (for mesons) in terms of contributions of 
higher Fock states. An appropriate framework for implementing these constraints was de- 
veloped in Ref. J7j: it is based on the derivation of EOM relations for non-local light-ray 
operators jTHj, which are solved order by order in the conformal expansion; see Ref. for 
a review and further references. In this way it is possible to construct self-consistent approx- 
imations for the DAs, which involve a minimum number of hadronic parameters. Another 
approach, based on the study of renormalons, was suggested for twist 4 in Refs. [2II1 EI] : this 
technique is appealing as it allows one to obtain an estimate of high-order contributions to 
the conformal expansion which are usually omitted. In this paper, we generalize this ap- 
proach to include S'f/(3)-breaking corrections and show how to combine renormalon-based 
estimates of "genuine" twist-4 effects with meson mass corrections. 

Pion DAs of twist 3 and 4 have already been studied in Ref. In Ref. [22]! these results 
were extended to the pseudoscalar octet; they include those meson-mass corrections that 
break chiral symmetry, while still preserving G-parity. S'[/(3)-breaking in the normalization 
of the twist-4 DA was estimated in Ref. ^3] • In this paper we continue the analysis of twist-3 
and 4 DAs of the K meson and present, for the first time, the complete set of DAs, including 
also G-parity-breaking terms that vanish in the limit of quarks with equal mass. The results 
are of direct relevance to the discussion of, for instance, 5-meson decays into light mesons 
using light-cone sum rules and also in the SCET framework. We refrain from an analysis of 
the 7] DAs, as the inclusion of the singlet part is crucial for phenomenological applications, 
but goes beyond the scope of this paper, and in fact deserves a separate study. 

Our paper is organized as follows: in Section |21 we explain notation and review the existing 
information on leading-twist DAs. Section |21 is devoted to twist-3 DAs: we give their general 
classification, calculate meson-mass corrections, perform a conformal expansion and formulate 
models in terms of a few non-perturbative parameters. A similar programme is carried out for 
twist-4 DAs in Section^ In Section El we present numerical results for all DAs and conclude 
in Section [Ul with a short summary and outlook. The appendices contain a collection of 
relevant formulas, in particular the QCD sum rules for the relevant twist-2, -3 and -4 matrix 
elements. 
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2 General Framework and Twist- 2 DAs 



Light-cone meson DAs are defined in terms of matrix elements of non-local light-ray operators 
stretched along a certain light-like direction z^, = 0, and sandwiched between the vacuum 
and the meson state. We adopt the generic notation 

0i;A/(M), A;m{u), ■■■ (2.1) 

and 

^t;M{a), ^t;M(«), ... (2.2) 

for two-particle and three-particle DAs, respectively. The first subscript t = 2,3,4 stands 
for the twist; the second one, M = it, K, . . ., specifies the meson. For definiteness, we will 
write most of the expressions for K mesons, i.e. sq bound states with q = u,d. The variable 
u in the definition of two-particle DAs always refers to the momentum fraction carried by 
the quark, u = Ug] u = 1 — u = Ug is the antiquark momentum fraction. The set of variables 
in the three-particle DAs, a = {01,0:2,03} = {o<j, Og, o^,}, corresponds to the momentum 
fractions carried by the quark, antiquark and gluon, respectively. 

To facilitate the light-cone expansion, it is convenient to introduce a second light-like 
vector such that 

P, = P,- (2.3) 
2 pz 

where is the meson momentum, = m\.j. We also need the projector onto the directions 
orthogonal to p and z, 

atu = 9t^u - —(Pi^^^ + P^^i^)^ (2-4) 

and use the notation 

az = a^z'', ap = af,p^, h^z = V^''' etc. (2.5) 

for arbitrary Lorentz vectors and tensors h^^y. a± denotes the generic component of 
orthogonal to z and p. 

We use the standard Bjorken-Drell convention [21] for the metric and the Dirac matrices; 



in particular, 75 = 'i7°7^7^7^, and the Levi-Civita tensor e^^x^ is defined as the totally 
antisymmetric tensor with €0123 = 1. The covariant derivative is defined as = — igA^ 
and the dual gluon-field-strength tensor as G^i, = ^e^^pa^C". 
The leading twist-2 DA of the K meson is defined as^ 

{0\q{z)[z,-z]^zl5s{-z)\KiP)) = tfKipz) ciw e^«^^ 02;x(u, /i') . (2.6) 



Here [x, y] stands for the path-ordered gauge factor along the straight line connecting the 
points X and y: 



[x, y]=P exp 



(2.7) 



ig dt{x-y)^A^{tx + {l-t)y) 
. Jo 

and n is the renormalization (factorization) scale. We also use the short-hand notation 

^ = 2m - 1. (2.8) 

^The leading-twist DA oi a, K meson is given by <p2-Ki'^) — 4>2:k{^ — u). 
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The decay constant fx is defined, as usual, as 



{0\mi,l5smK{P)) = iIkP,, (2.9) 

with fx = 160 MeV jHl- The normahzation follows from the requirement that the local limit 
z of reproduce (jTHl), so that 

! du(P2;K{u) = 1. (2.10) 

Jo 

A convenient tool to study DAs is provided by conformal expansion JH].^ The underlying 
idea is similar to partial-wave decomposition in quantum mechanics and allows one to separate 
transverse and longitudinal variables in the Bethe-Salpeter wave-function. The dependence 
on transverse coordinates is formulated as scale dependence of the relevant operators and is 
governed by renormalization-group equations, the dependence on the longitudinal momentum 
fractions is described in terms of irreducible representations of the corresponding symmetry 
group, the coUinear conformal group SL(2,R). The conformal partial-wave expansion is 
explicitly consistent with the equations of motion since the latter are not renormalized. It 
thus makes maximum use of the symmetry of the theory to simplify the dynamics. 

To construct the conformal expansion for an arbitrary multiparticle distribution, one first 
has to decompose each constituent field into components with fixed Lorentz-spin projection 
onto the light-cone. Each such component has conformal spin 

J = l{^ + s), 

where / is the canonical dimension and s the (Lorentz-) spin projection. In particular, 1 = 3/2 
for quarks and / = 2 for gluons. The quark field is decomposed as ip^ = A^ip and ip^ = A_il) 
with spin projection operators A+ =^^7 (2pz) and A_ = ^/ {2pz), corresponding to s = +1/2 
and s = —1/2, respectively. For the gluon field strength there are three possibilities: Gz± 
corresponds to s = -|-1, Gp± to s = —1, and both G±± and Gzp correspond to s = 0. 
Multiparticle states built of fields with definite Lorentz-spin projection can be expanded in 
irreducible representations of SL(2,]R) with increasing conformal spin. The explicit expression 
for the DA of an m-particle state with the lowest possible conformal spin j = ji + ■ ■ ■ + jm, 
the so-called asymptotic DA, is p/7| 

0a.(ai, "2, ■ ■ ■ , am) = ^^^7p-V — FT^T^ al^'~^al^'~^ . . . 0^""^ (2.11) 

r(2ji) • ■■T{2jm) 

Multiparticle irreducible representations with higher spin j + n,n = 1,2, . . ., are given by 
polynomials of m variables (with the constraint X^fcLi Q^fc = 1 ); which are orthogonal over 
the weight function (|2.11|) . 

In particular, for the leading- twist DA (f)K;2 defined in fl2.6p . the expansion goes in Gegen- 
bauer polynomials: 

0X;2(m, /i') = 6uil -u)(^l + f2 a^ifi^)Cl/\2u - 1) j . (2.12) 
^See Ref. |25| for an alternative approach not based on conformal expansion. 



4 



To leading- logarithmic accuracy (LO), the (non-perturbative) Gegenbauer moments a„ renor- 
malize multiplicatively with 



-7i°V(2/3o) 



«n(/^o) 



(2.13) 



where L = as{fi'^)/as{f4i), 
by 



[llNc — 2Nf)/3, and the anomalous dimensions "fn^ are given 



8Cf 



(2.14) 



2(n + l)(n + 2), 

The reason why leading-order renormalization respects the (anomalous) conformal symmetry 
is that it is driven by tree-level counterterms that retain the symmetry properties of the La- 
grangian. More technically, the Callan-Symanzik equation that governs the renormalization- 
scale dependence coincides to this accuracy with the Ward identity for the dilatation operator, 
which is an element of the coUinear conformal group [T^ . 

To next-to-leading order (NLO) accuracy, the scale dependence of the Gegenbauer mo- 
ments is more complicated and reads j2ni 



,NLO 



k=0 



where 



7^^'' are the diagonal two- loop anomalous dimensions |23, (3i = 102 — (38/3) A^/, and the 



E, 



NLO 



^7i°V(2/3o) 



1 + 



(2.15) 



(2.16) 



mixing coefficients S^^^ k < n — 2, are given in closed form in Ref. i27j, see also, for instance. 



Ref. 



for a recent compilation. For the lowest moments = 0, 1, 2 one needs 



76 



(1) 



7i 



(1) 



23096 512 



243 



N 



f ' 



34450 830 



243 



Nf 



and 



(1) 

20 



/5o) 



(0) 
72 



JO) 



2^-1+7^ /(2A,)' 



(2.17) 



(2.18) 



72^' - Wo 

The odd Gegenbauer moments a2„+i are first order in S'?7(3)-symmetry breaking for the 
kaon and vanish for the pion by virtue of G-parity. The numerical value of was the 
subject of significant controversy until recently. The existing estimates are all obtained using 
QCD sum rules. The first calculation of by Chernyak and Zhitnitsky yielded af ~ 0.1 
Pm 1^ . but unfortunately suffers from a sign mistake in the perturbative contribution [7j. 
After the error is corrected, one finds that the two numerically leading contributions come 
with different sign and cancel to a large extent, so that the sum rule becomes unstable and 
numerically unreliable. This problem was reanalysed in Refs. jHl IHl CSl using a different 
set of sum rules, where it was also checked that the results are consistent with the equations 
of motion for the relevant operators [HI E] • The results are given in Table ^ As our best 
estimate, we take 



IGeV) = 0.06 ±0.03. 



(2.19) 
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Method 


^ = 1 GeV 


/i = 2 GeV 


Reference 


QCDSR,D 


0.05 ±0.02 


0.04 ±0.02 


i 


QCDSR,ND;EOM 


0.10 ±0.12 


0.08 ±0.09 


i 


QCDSR,D;EOM 


0.06 ±0.03 


0.05 ±0.02 


dmn] 



Table 1: The Gegenbauer moment af^(/i^) from QCD sum rules. The abbreviations stand for: 
QCDSR: QCD sum rules; D and ND: diagonal and non-diagonal correlation function, respectively; 
EOM: equations of motion. The error estimates should be taken with some caution, as there is 
no systematic approach to estimate uncalculated higher-order terms in the OPE. 



Method 


fi = 1 GeV 


fi = 2 GeV 


Reference 


CZ model 


0.56 


0.38 




QCDSR 




1 7+0- 14 

-0.06 


E 


QCDSR 


0.28 ±0.08 


0.19 ±0.05 


this paper 


QCDSR,NLC 


0.19 ±0.06 


0.13 ±0.04 


jsaisaiHH] 




0.19 ±0.05 


0.12 ±0.03 (/i = 2.4) 


m 




0.32 


0.20 ifx = 2.4) 


m\ 


Fj^jj* ,LCSR,,R, 


0.44 


0.30 






0.27 


0.18 


m 


f;"^,lcsr 


0.24 ±0.14 ±0.08 


0.16 ±0.09 ±0.05 


gnum 


F;^"^,LCSR,R 


0.20 ±0.03 


0.13 ±0.02 


m 


-^B-^TT^i/^LCSR 


0.19 ±0.19 


0.13 ±0.13 


m 


LQCD, quenched, 


0.381 ± 0.2341°:^^^ 


0.233 ±0.143l[]:|]| 


UKQCD 


w/cw 




ifi = 2.67) 




LQCD, Nf = 2, W/CW 


0.364 ±0.126 


0.236 ± 0.082 = 5) 


QCDSF/UKQCD fill 



Table 2: The Gegenbauer moment al^^i^). The CZ model involves al^ = 2/3 at the low scale 
fx = 500 MeV; for the discussion of the extrapolation to higher scales, see Ref. [37j. The 
abbreviations stand for: QCDSR: QCD sum rules; NLC: non-local condensates; LCSR: light-cone 
sum rules; R: renormalon model for twist-4 corrections; LQCD: lattice calculation; Nf = 2: 
calculation using Nj = 2 dynamical quarks; W/WC: Wilson glue and non-perturbatively 0{a) 
improved Clover-Wilson fermion action. 
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Calculations of the second Gegenbauer moment for the pion DA, Og, have attracted quite 
a bit of attention and have a long history. Three different approaches have been used: direct 
calculations using QCD sum rules, pioneered by Chernyak and Zhitnitsky; analysis of the 
experimental data on the pion electromagnetic and transition form factors and the B weak 
decay form factor, using light-cone sum rules; and lattice calculations. The summary of 
these results is presented in Table 12 see also, for instance, Refs. |37| 1^ for another recent 
compilation. 

Our conclusion from Table |21 is rather pessimistic: can only be determined with large 
errors, whatever approach is chosen. A fair quote is probably 

4(1 GeV) = 0.25 ±0.15. (2.20) 

The i^-meson DA has attracted comparatively less attention. The old estimate by 
Chernyak and Zhitnitsky, {i"^) k / {i"^) ^ = 0.8 ± 0.02 ^Tj, translates to 

= 0.59 ±0.04 ^ (af )cz(l GeV) = 0.33 (2.21) 

for the CZ model. A recent calculation, Ref. [H], including radiative corrections to the sum 
rules gives, however 

af /4 ~ 1 , (1 GeV) = 0.27t'oil. (2.22) 

This result is consistent with the most recent lattice calculation, using Nj = 2 dynamical 
fermions [44,; which shows that (^^)7r stays practically constant under a variation of the pion 
mass. For the purpose of the present paper we have done an update of the QCD sum-rule 
calculation |H5, using the corrected 0{as) quark-condensate contribution given in Ref. |10j . 
see App. |BJ which yields 

af/a^ = 1.05 ±0.15, af(l GeV) = 0.30 ± 0.15 . (2.23) 

The difference with \^ is small and mainly due to the larger value of the strong coupling 
that we are using in this work. We conclude that the existing evidence points towards a 
very small SU{3) violation in the second coefficient in the Gegenbauer expansion, so that we 
accept af' = in the range given in Eq. (j2.20p as our final estimate. 

Estimates of yet higher-order Gegenbauer coefficients are rather uncertain. The light-cone 
sum-rule calculations of the transition form factor F^r^^. in Refs. ^Hl EH IHHl (HH] suggest a 
negative value for al, which is consistent with the result aJ(lGeV) > —0.07 obtained in 
Ref. pni- However, this conclusion may be premature because of the omission of yet higher 
moments and absence of any convincing method to estimate systematic errors involved in the 
analysis. For this reason we adopt, in this paper, a model for the leading-twist DA, which is 
given by the Gegenbauer expansion ()2.12p truncated after the second term. 

Last but not least we have to specify the value of the strange-quark mass. In the last year 
several lattice calculations with dynamical fermions have been published; see Refs. |45, 46] 
for a summary and short review. In all these calculations the physical kaon mass is used as 
an input to fix the strange-quark mass. There is good agreement between data sets obtained 
using different non-perturbative renormalization procedures and, in fact, also with earlier 
quenched calculations. On the other hand, the data still show considerable dependence on 
the lattice spacing, so that it is clear that simulations on finer lattices are needed for a 
systematic continuum extrapolation. In a different approach, the strange-quark mass can be 
extracted from the e~^e~ annihilation cross section and/or hadronic r-decay data using QCD 
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sum rules. These calculations have reached a certain degree of maturity and yield results that 
are in reasonable agreement with lattice determinations; see Ref. ^| for a recent summary 
and further references. In this paper we use 

m,(2 GeV) = (100 ± 20) MeV, (2.24) 

which corresponds to = (137 ± 27) MeV at 1 GeV. 



3 Twist-3 Distributions 

In this section we define all the twist-3 DAs of the kaon and derive models for them to 
next-to-leading order in conformal expansion, which fulfil the QCD equations of motion. We 
also work out the leading-order scale-dependence of the corresponding hadronic parameters. 
Numerical values for the parameters and the corresponding models are given in Section El 
To twist-3 accuracy, there are two two-particle DAs defined as 

m{^)^l^s{-z)\K{P)) = -^^^ /'rfne^(2"-i)^^0f.^(«), (3.1) 

ms + mq Jo 

{0\q{z)a^p^,s{-z)\KiP)) = J^Sl^{P^zp - PpZ^) f du^^'^-'^^^ r.-A^)- (3.2) 

In addition, there is also one three-particle DA: 

{0\q{z)a^,^59Ga(s{vz)s{-z)\K{P)) = 

= I fsK {PaP.g^p - PaPug^p - (« ^ /?)) j P« e-^^^("^-°^+^"^) $3;i^ («! , ^2, ^s) + • 



where the integration measure is defined as 



J T>a = J daida2das5{l — ai — a2 — ce^) 



(3.3) 



(3.4) 



and the dots stand for Lorentz structures of twist 5 and higher. 
To next-to-leading order in conformal spin, $3;;^ is given by 

^3;Kia) = 360aia2al |l + A3ii-(ai - 02) + ujsk ^ (7^3 - 3)| . (3.5) 

The three parameters f^x, ^sk, and u^k can be defined in terms of matrix elements of local 
twist-3 operators as follows: 

{0\qa^(.-f5gG^i.s\K) = 2ifsKipzf, 

3 3 

{0\qcr^i.'j5[iD^,gG^^]s - -id^qa^^.'j^gG^i.slK) = 2if^K{pzf — ^kz^ 

{OlqiDzCrz^^gGz^s - qa ^^-i^^gG z(,i D z s\K) = 2if^K{.pzf ^ Xk'a- (3.6) 
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Numerical values for these parameters can be obtained from QCD sum rules and will be 
discussed in Section |S] 

The operators in ()3.6|) renormalize multiplicatively in the chiral limit, with one-loop 
anomalous dimensions jlSj 

73;/ - ^"-^^ Y ^ 

ro) _ 20 7 208 



0:i{z,vz,^Y = 02{z,vz,QY° -ims-^\\i^- / dt 02{z,vz) - 2t02{z,tvz) , (3.8) 



3 " ' 3^ 9 ' 

,<«..ic. + fc.^H?. (3.7, 

For a massive strange quark, the operators in ()3.6|) mix with twist-2 ones. Using the light- 
ray-operator technique of Ref . ^H] 5 this mixing can be expressed in compact form as 

2tt V Jq 

where 

Os{z, vz, of = [qiz)a,,^59G,,ivz)s{0)f (3.9) 

and 

02(0^,^)^^' = [q{azY,^,s{bz)f; (3.10) 

/i^ stands for the normalization point. Sandwiching (|3.8|) between the K state and the vac- 
uum, and expanding in powers of p^, one can easily derive the mixing for local operators with 
an arbitrary number of derivatives. We find that fsx mixes with fx^T^s and with fxragai , 
whereas X^k and u^k mix in addition with fxmsa2- The corresponding LO renormalization- 
group-improved expressions read 

+ ^ (^55/(9/30) _ ^68/(9/30)) ^^^^^^i.](^2)^ 
^ ^ (^68/(9/3o) _ ^139/(18A0) [rri.af ] (/x^) 
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where L is the leading-log scaling factor: L = Q;s(/i^)/as(yWo)- 

The two-particle twist-3 DAs ()3.1|) and ()3.2|) are not independent, but related to the 
three-particle DA $3;A' by EOM ^^lESj- The EOM relations contain terms that depend on 
quark masses and can conveniently be expressed in terms of two dimensionless parameters 

^+ " ml ' ml ' ^^-^^^ 

numerically ^ p^. The rationale for introducing two parameters is that p5 changes sign 
when switching from K mesons to K mesons, i.e. p^ = p^, but p5 = —P-- In the analysis 
of twist-3 DAs given in Ref. j22]; only terms in p^ have been included. Here we complete 
these studies by taking into account also the terms in p^. 

From the non-local operator identities ()A.3|) and ()A.4|) . one obtains the following relations 
for moments of the DAs, dropping the index K: 

n + 6 n + 6 n + 6 

-P+^ Mt +P-^ Mt, , (3.13) 
where we use the notation 

r-l 



M^= du{2u-l)"(l){u) 
Jo 

and introduce the auxiliary functions 

I'll ^2 

ipil\u)= dai da2- $3(«), (3.14) 

Jo Jo l-ai- a2 

ipf\u)= dai da2- («i - «2 - (2« - 1)) $3(«) . (3.15) 

Jo Jo {l-ai-a2y 

The normalization is chosen in such a way that 

1 1 

= I du(pl{u) = 1, M^'' = [ ducp^iu) = l-p+. (3.16) 



Except for the new terms in p_, these moment relations agree with those obtained in Refs. ^3 
123. 

The relations ()3.13p can be solved exactly: separating the contributions of quark-anti- 
quark-gluon operators and the terms in p±, 

4>l{u) = 1 + <Pl^{u) + p+ <j>l4u) + p^ <Pl_{u) , 

r^iu) = 6uu + r3,g{u) + P+ 03,+(«) + 03,-(«) , (3.17) 

10 



we obtain the integral representations (cf. Ref. 



1 rdv 



4 ./n V 



,(2)^ 



1 r^di; 



4 / f 



{2v - 1) (v.^^))"(t;) - 2 (v.rO'(^) + (V^r)"(^) , (3.18) 



,(1) 



(2) a"/ 



03,+ («) 



4io 4J„ f 



(3.19) 



J r 1^ [202(t;) - Uv)] - I ^ + <^2(^)] , (3.20) 

^ Jo ^ ^ Ju ^ 



where primes denote the derivatives in v. (p'iy) = {d/dv)(l){y) etc. 

For the second twist-3 DA the solutions of the moment relations read, in the same nota- 
tion: 



1 2 



V 



KM = -^uuU dv(- + ^\{ (^«)'(^) + {2v - 1) {vf^)'{v) 



(3.21) 
(3.22) 



^3,+ (^) = dv^(f)2{v) + j dv^(f)2{v)^ 



'^Ju) = luulf dv f 4 + ^] Mv) - ! rff f 4 + - I Mv] 



2 1 /n Kv"^ V L, \ f 2 f 



(3.23) 



We stress that the relations ()3.18j) to ()3.23j) are valid in full QCD and involve no approxi- 
mation whatsoever. One consequence of these relations is that quark-mass corrections to 
contain logarithmic end-point singularities. In particular for the asymptotic leading-twist DA 
(p2;K{u) = 6u{l — u) we obtain 



0l;i^(«)| 
Ir 



1 + p5 ^ (2 + Iiimm) + p:^ - f 1 - 2m + In -V (3.24) 
2 2 V u/ 



= 1 + 

I no gluons, asymptotic f/>2;X 

To NLO in conformal spin we obtain, using the truncated conformal expansions ()2.12|) for 
(t)2-K and (|3.5p for ^z\K'- 



.j,{u) = 1 + 3p^(l + 6af ) - 9p^af + Cl'\2u - 1) 



+Cl'\2u - 1) (30r/3i^ + ISp^af - 3p5af ) + Cl'\2u - 1) {iQ^k^^k - \ p'^af^ 



3r]3KUJ3KCy\2u - 1) + - (p^ + p^)(l - 3af + 6af ) \nu 
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+ 2 (Z^? - P"") (1 + 3«f + 64)111 



(3.25) 



'u) = 6uu 
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+ {^VzK - \ mKU^3K + I P+af ^ C|/'(2u - 1) + r]sKX3KCl^\2u - 1) 



+ 9uu{pI +p^){l-3af +Qa!^)\nu + 9uu{pl -p^){l + 3af +6a^)ln 
where, to simplify notations, we have introduced the parameter 



u . 



V3K 



Ik m 



(3.26) 



(3.27) 



K 



These expressions are our final results for the two-particle twist-3 DAs and supersede those 
given in Refs. ^Tj The terms multiplying Inw and Inw are the first three terms in the 
conformal expansion of 02 i<-(O) and 02.^(1), respectively. Numerical values for the hadronic 
parameters are given in Table El The leading-order scale-dependence follows from 1)3.111) and 
the scale dependence of the quark masses in p^ and rj^K- 
We note in passing that the EOM relations 



no gluons 



no ffluons 



no ffiuons 



K 



are no longer fulfilled for (p^.^, but violated by mass corrections in p^ 



4 Twist-4 Distributions 

In this section we derive models for the two- and three-particle twist-4 DAs to NLO in 
the conformal expansion. There are four i^-meson three-particle DAs of twist 4, defined 
as [1711221^ 

m{z)l,1^9Gc.p{vz)s{-z)\K{P)) = 

= Pf^iPaZp - PpZa) — fK^i;K{.V,pz) + {ppQ^^ - p^g^^) f ^-k^V ,pz) + . . . , (4.1) 

pz 

{mz)l,wGo.p{vz)s{-z)\K{P)) = 

= P^i{PaZp - PpZa) — fK^4;K{v,pz) + {ppQ^^ - Pag^^,) f i;K{v , pz) + . . . , (4.2) 

P£ 

^In the notation of Ref. [221, ^i;K = ™/f"^lh *4;^ ^ 'mj^A±, $4;A' = and *4;/f = mj^V±. 
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with the short-hand notation 



The integration measure Va is defined in (j3.4|) . and the dots denote terms of twist 5 and 
higher. For massless quarks and, more generally, for two equal-mass quarks, G-parity implies 
that the DAs $ and \l/ are antisymmetric under the interchange of the quark momenta, 
ai ^ a2, whereas $ and \l/ are symmetric pT| l22j. Note that unlike twist-2 and twist-3 DAs, 
which are dimensionless, the twist-4 DAs have mass dimension 2 (GeV^). The corresponding 
contributions to hard exclusive processes are suppressed by two powers of the hard scale with 
respect to leading twist. 

The distribution amplitudes and correspond to the light-cone projection ■jzGzp, 
which picks up the s = +1/2 components of both quark and antiquark field and the s = 
component of the gluon field. The conformal expansion reads: 



$4;i^(a) = 12Oaia2a3[0^ + 0f ("i - "2) + 0f (Sas - 1) + • • •] 
$4;i^(a) = 12Oaia2a3[0^ + 0f ("i - "2) + (Saa - 1) + • • •] 



(4.3) 



G-parity implies that, for the vr meson, 0q = 02 = 0i = O5 whereas 0^, 0^ and 0f are 
0{ms-mq). 

In turn, the DAs ^^-^k and '^^-k correspond to the hght-cone projection 7_lG'^_l, which is a 
mixture of different quark-spin states with Sq = +1/2, Sq = —1/2 and Sg = — 1/2, Sq = +1/2, 
respectively. In both cases s = -|-1 for the gluon. We separate the different quark-spin 
projections by introducing the auxiliary amplitudes and defined as 



The original distributions "^a-k and "^^-k are given by 



2 - 



(4.4) 



(4.5) 



"^^^ and have a regular expansion in terms of conformal polynomials, to wit: 



vl>Ti| 



a 



600203 



i'o^ + ^{^("3 - 3ai) + V^2^ ( as - 7^a2 



Ti 



i'o + - 3^2) + ^2' ( «3 - ^ttl 



it 



For the tt meson, thanks to G-parity, 
so that ipj^ = '01^.'^ For K, we write 

^j^ = i,f + ef, = - ef. 



(4.6) 



(4.7) 



(4.8) 



*This implies, in particular, that only one of the DAs '5 and ^' is dynamically independent. 
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where the 6i correspond to S'f/(3)-breaking corrections that also violate G-parity. From ()4.5|) . 
the following representations can readily be derived: 



^4;7^(«) = -30q;3<{ V'o' (1 - "3) + V'f "3(1 - "3) - 6a;ia2 



K 



(4.9) 



03(1 - as) - -{al + al) 



- (ai - ^2 



^0 + «3^f + - (5a3 - 3)^2 



K 



^4;i^(a) = 30a2|^^(i - 03) + 03(1 - 03) - 6aia2 



(4.10) 



K 



a^{l - a^) - -{al + al) - (ai - ^2) 



1 



+ as^Af + - (5^3 - 3)^: 



K 



In addition, we introduce one more three-particle DA S4(«) |21j : 



(0 |g(^)7^75^7l^°G„/3(t^^)s(-^)| K+(p)> = ^/xp^p^ y P« e-*^ 



ipz{»2-a,+v^^) ^^^^^^^y (4.11) 



The Lorentz structure p^pp is the only one relevant at twist 4. Because of the EOM, D^G^^ = 
~9Yl,q9'^lp^^ where the summation goes over all light flavors, S4;ft:(a) can be viewed as 
describing either a quark-ant iquark-gluon or a specific four-quark component of the pion, 
with the quark-ant iquark pair in a color-octet state and at the same space-time point. The 
conformal expansion of 24.^^(0) starts with J = 4 and reads 



840aia2«3 



-t- . . . 



(4.12) 



where has mass dimension 2. The dots stand for terms with higher conformal spin 
J = 5, 6, . . ., which are beyond our accuracy. This DA was not considered in Refs. [T71 

= 0{ms — rriq) and vanishes for mesons built of quark and antiquark with equal 



because 
mass. 

The expressions in Eqs. (14. 3|) . (j4.9p . (j4.1(Jj) . (j4.1ip represent the most general parametriza- 
tion of the twist-4 DAs to NLO in the conformal-spin expansion and involve 13 non-perturba- 
tive parameters. Not all of them are independent, though. In the following, we shall establish 
their mutual relations and also express the expansion coefficients in terms of matrix elements 
of local operators. 

The asymptotic three-particle DAs correspond to contributions of the lowest conformal 
spin J = js + iq + = 3. The parameters 0^, 0^, ip^ and 6^ describing these DAs can be 
expressed in terms of local matrix elements as 

mic.gG^c.s\K{P)) = iP^fK^l, 



(4.13) 



These are the only two local twist-4 operators of dimension 5. Note that the second matrix 
element vanishes for equal-mass quarks, because of G-parity. It also vanishes in the chiral 

Moreover, in this limit can be calculated 



limit mq, 



because of the factor m\. 



exactly to leading order in 0: 



(4.14) 
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numerical estimates of the corrections can be obtained from QCD sum rules and will be 
discussed below. 

Taking the local limit of Eqs. ()4.H1 . ()4.2|1 . one obtains 

< = ^o=-l 4, <l>o = = I ^Wk- (4.15) 

What about the scale-dependence of these parameters? Like f^x, renormalizes multi- 
plicatively for massless quarks, but mixes with operators of lower twist for 7^ 0. At the 
operator level, neglecting 0{m1) corrections, the mixing is given by 

{qiagG^o^sf = {qjagG,.sf^ ("i _ ^ ^ In 4 V ^ - In 4 mMT' ■ (4-16) 

Taking matrix elements and resumming the logarithms, we find 

[6],m = + ^ (1 - L^'^'''^') mj,, (4.17) 

with, as before, L = as(/i^)/a;s(/io). 

The scale dependence of k,4k can most easily be derived by observing that this parameter 
is related to af and quark masses by the equations of motion 0: 

1 iris - rria 9 ml - ml , s 

K^^ = f i af+ ^ ^ 4.18 

Taking into account the known scale dependence of af" and g, one obtains 

^Mf^') = ^.M) - I {L'^'^''^' - 1) af (/.^) + [L^'^^ - 1) ^"^-^^^^^^"^ . (4.19) 

To NLO in conformal spin, the discussion becomes more involved. As explained in 
Ref. P7j, for massless quarks the corresponding contributions can be expressed in terms 
of matrix elements of the three existing G-parity-even local quark-antiquark-gluon opera- 
tors of twist-4. These three operators are not independent, however, but related by the QCD 
equations of motion. One is left with one new non-perturbative parameter only, call it uj^k,^ 
which can be defined as 

4 _ _ 

(0|g[iL'^, igGy^]-i^s - - id^qigGy£_-i^s\K{P)) = 

= fK^loJiK (^P^Pu - \ mlg^^j + 0(twist 5). (4.20) 
The scale dependence of uj^k, for massless quarks, is given by 

For massive quarks, a distinction must be made between G-parity-conserving and G-parity- 
breaking contributions. G-parity-conserving corrections do not involve new operators, and 



^In the notation of Ref. ^0Ji = (8/21)e. 
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the difference to the massless case is mainly due to corrections proportional to the meson 
mass. This case is described in detail in Refs. j2nil221- Here we just quote the results obtained 
in Ref. 123: 



,2 „K 



7 ,2 3 



21 2 

^2 = ^K^iK, 



3 

20 



K^2 ■ 



(4.21) 



The G-parity-breaking contributions, on the other hand, involve a different set of local oper- 
ators and in particular 



which determines the normalization and the leading conformal spin contribution to the DA 
S4;i^(a) defined in Eq. ()4.11|) . Hence, a complete treatment of G-parity-breaking corrections 
to twist-4 DAs requires also the inclusion of 'B^-k- 

It is beyond the scope of this paper to work out the corresponding relations between 
the matrix elements of local operators and expansion coefficients. For this reason, and also 
because QCD sum-rule estimates of matrix elements of large mass dimension are not very 
reliable, we adopt a different approach and estimate G-parity-breaking corrections of spin J = 
4 using the renormalon model of Ref. [21] . The general idea of this technique is to estimate 
matrix elements of "genuine" twist-4 operators by the quadratically divergent contributions 
that appear when the matrix elements are defined using a hard UV cut-off, see Ref. |^ for 
details and further references. In this way, three-particle twist-4 DAs can be expressed in 
terms of the leading-twist DA 02 : 



^4':^ (01,02,03) 



r)2 



1 — Oi 



h;K[a2) 
1 - «2 



$^^^(01,02, as) 



a2<P2;K[ai, 

(l-«l)2 



<yi<P2;K{<y2) 
(l-«2)2 



^'4':^(ai,a2,a3) 
$4^^(ai,a2,a3) 



S^^^(ai, 02,03) 



6 



261 



(f)2;K{ai) _^ (j)2;K{a2) 



1 — ttl 1 — 0:2 



a2<p2;K[ai) _^ ai02;i^l«2j 



3 L (l-«i)' (l-«2)' J 

a2 02;if(ai) ai02;/^(tt2) 



1 — «! 



1-02 



(4.22) 



where, in difference to [21], we do not assume that 02 (w) is symmetric under the exchange 
u ■(-^ 1 — u. 

The expressions in (j4.22j) do not rely on conformal expansion and contain the contributions 
of all conformal partial waves. Projecting onto the contributions with the lowest spin J = 3, 4, 



16 



we obtain 



C = 0, e^< = ^af6l, e^ = -\4b\. (4.23) 
It follows that in the renormalon model 

2 

U)^K = t^47r = -, (4.24) 

y 

which is in good agreement with direct QCD sum-rule calculations ^7]. We also find 

H^ = iaf4. (4.25) 

Note that in the renormalon model = 0. This is due to the fact that the contribution 
in in Eq- ()4.15p is obtained as the matrix element of the operator ()4.13|) which vanishes 
by the EOM (up to a total derivative), see Eq. ()4.18|) . Therefore, against appearances, 
this contribution has to be interpreted as "kinematic" power correction induced by the non- 
vanishing i^'- meson mass rather than a "genuine" twist-4 effect. 

We are now in the position to derive expressions for the two-particle DAs of twist 4. They 
are defined as 

(0|g(x)[x,-x]7^75s(-x)|K(P)) = i/i^P^ [ due'^^'' i<p2Au) + \ x^(l)i,K{u) 







+ lf^J^^f^l' e^«^"^4;i.(w), (4.26) 



which is the extension of Eq. ()2.6|) to twist-4 accuracy.^ From the operator relations ()A.1|1 
and ()A.2|) . we obtain 



i^K;AW = m^|205.^(M) - 02;A'(m)| + 3" / "«1 / ««2 , (4.27) 

du Jq Jo 1 - ai - a2 

= 12^4;k{u) - 12m^</)2;x(n) - 2^ [{2u - l)(m|02;x(M) + ^4;k{^))] 



r, f, 4(2vl>4;K(«)-$4;i.(«)) , ... 

+— / dai / da2 rf (ai - - {2u - 1)) 

du^ Jo Jo (1 - ai - 02)^ 



ms - ma 9 d(jfn.j^{u) 
+ A— -mi 4.28 

ITLs + mq du 



^fpA;K and (j)4-K are related to the DAs defined in Ref. by 04;/f = m\gK and ?/'4;i<: = m\^KK- 
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with the boundary condition 04;ic(O) = 04;k(1) = 0. 

We solve this relations splitting the result in "genuine" twist-4 contributions ipj.j^ and 
Wandzura-Wilczek-type mass corrections ip^jf as 

^4;i.(w) = ^IMu) + (4.29) 

with 

i^l'M = Y ^kCTK-^u - 1) + 5 {Sef - } Cy\2u - 1) , (4.30) 
^J7(„) = {1 + 6p^(l + 6af ) - ISp^af } Co/'(2m - 1) 

+ m\ |-12/€4i^ - ^ «f + 27p:^af - 3p:^(l + 18af )| Cy^(2?i - 1) 

+ (^1 + y af + aOp^af - 6p5af ^ + 60 ^ (m, + m,)| C2'/'(2n - 1) 

+ {< + y '^4/. - 9P~4) + 20 ^ (m. + m,)A3i^| C]!\2u - 1) 

+ |~^"^^f'^2^ 6 ("^^ + "^'?)^3i^| C'4^^(2m - 1) 

+ Qmg{ms + 'mg){l + Saf' + 6a^) In-u + 6ms(ms + rng){l — 3af" + 6a^) Inu , 



(4.31) 



where ^ = 2u — 1, see Eq. ()2.8p . ip^x" vanishes for rriK — * and ms,q — > 0. 
The complete expression for 



M^) = <Pl%{n) + <P^^{u) (4.32) 
is rather lengthy. We find for the "genuine" twist-4 part: 

+ 216]^uj4k[uu{2 + 13mm) + [2m^(6m2 - 15m + 10) Ium] + [m m]| 

+ 4O0f {mm^(2 - 3mm) - [2m^(m - 2) Ium] + [m ^ m] } , (4.33) 
and for the mass-corrections, neglecting numerically small terms of order m^: 

m|^K4;^|MM^(l — 2mm) + [5(m — 2)M'^lnM] — [m ^ m]| 
+ 4^(m, + mg)MM<|30 ( 1 ' ^ 



iWW, 



M 



fx { \ ms+TTlq 
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+ lOAsK e [1 



uu\ 



rris + rrio 



[1 — 5uu] 



[3- 21mm + 28u^u^ + 3^ 



m.o - rua 



rris + mg 



[1 - 7mm] 



ym|-af |^mm(4 - 9mm + IIOmV) + [m^(10 - 15m + 6m^) Ihm] + [m ^ m]| 



+ 4m% uu (1 + 3mm) ( 1 + ^af ^ ) . 



(4.34) 



The DAs for K mesons are obtained by replacing m by 1 — m. Note that ip4,-^K has logarithmic 
end-point singularities for finite quark mass, whereas (p^.^K has no such singularities, so that 
one can safely neglect the 0{m1) terms. 

The expressions given above provide a self-consistent model of the twist-4 DAs which 
includes the first three terms of the conformal expansion.^ An estimate of the contribution of 
higher orders can be obtained using the renormalon model. In this case, the "genuine" twist-4 
contributions to the two-particle DAs given in Eqs. ()4.29j) and fl4.32j) have to be replaced by 



,r4,ren/ 



M 



-8 



K 



1 r 1 

dv(j)2-K{v) \ ^ 



M + M 



m) In 



1-^ 

V 



e{v-u) 



+ 



u + M + (m — m) In ( 1 — 



-\]eiu~v) 
v) J 



/T4,ren/ \ 



81 l_ 

3 du?' 



dv 4>2;k{v) 



e{v-u) + 

V / 



M^ ^ 



(4.35) 



and used in combination with the complete renormalon-model expressions for the three- 
particle DAs given in Eq. ()4.22p . As explained in Ref. the renormalon model does 
not take into account the damping of higher conformal-spin contributions by the increas- 
ing anomalous dimensions and, therefore, provides an upper bound for their contribution. 
The effect of these corrections is, most importantly, to significantly enhance the end-point 
behaviour of higher-twist DAs in some cases, which can be important in phenomenological 
applications. 



5 Models for Distribution Amplitudes 

In this section we compile the numerical estimates of all necessary parameters and present 
explicit models of the twist-3 and -4 two-particle distribution amplitudes that we introduced 
in Sections El and 01 The important point is that these DAs are related to three-particle 
ones by exact QCD equations of motion and have to be used together; this guarantees the 
consistency of the approximation. Our approximation thus introduces a minimum number of 
non-perturbative parameters, which are defined as matrix elements of certain local operators 

''One shortcoming of the model is that G-parity-breaking meson mass corrections of spin J = 4 are missing 
and we only include the "genuine" G-parity-breaking twist-4 corrections estimated in the renormalon model. 
Numerically, both effects may be of the same order. 
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taKen irom rtei. jiu] . vor-oQci 




U.ZO IC U. iO 
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u.zo m u. io 
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U. i ( IC u. lU 


/ / 1 -4l r^ivdol^in (T L ' y > "1 oil' 

0(7 1^0 J UIcclKlIlg bllidii, 

average over varioub Refb. 


J oi\ 


0.45 ± 0.15 


0.33 ± 0.11 


J OTT 


0.45 ± 0.15 


0.31 ± 0.10 


in unitb of 10'^ GeV^ 




-1.2 ±0.7 


-0.9 ±0.5 




-1.5 ±0.7 


-1.1 ±0.5 






1 « _)_ n 4 


1 4c; _)_ n 

X.HzU _1_ \J.OO 




n 

U 


n 


VJT tJU-U. 




0.20 ±0.06 


0.17 ±0.05 


51 


0.18 ±0.06 


0.14 ±0.05 


in unitb of GeV^ 




-0.09 ±0.02 


-0.10 ±0.02 










taken from Ref. G-odd 


UJ4K 


0.2 ±0.1 


0.13 ±0.07 




0.2 ±0.1 


0.13 ±0.07 


taken from Ref. |17j: 
SU (3) breaking not incl. 



Table 3: Hadronic parameters for the K DAs. We also give the corresponding parameters for 
the 71, which are a by-product of our calculations. All parameters have been calculated in this 
paper at n = 1 GeV, unless stated otherwise. The evolution between 1 and 2 GeV is done at NLO 
accuracy for frig s and a^'^ , and at LO accuracy for the other parameters. The twist-4 parameters 
6f, (pf etc. are given by Eq. ()4.23|) . based on the renormalon model. 

between the vacuum and the mebon btate, and which we ebtimate ubing QCD bum ruleb. 
More bophibticated modelb can be conbtructed in a bybtematic way by adding contributions 
of higher conformal partial waveb when ebtimateb of the relevant non-perturbative matrix 
elementb will become available. 

Our approach involveb the implicit abbumption that the conformal partial wave expanbion 
ib well convergent. Thib can be jubtified rigoroubly at large bcaleb, bince the anomalous 
dimensions of all involved operators increase logarithmically with the conformal spin J, but 
is non-trivial at relatively low scales of order /x ~ (1-2) GeV which we choose as reference 
scale. An upper bound for the contribution of higher partial waves can be obtained from the 
renormalon model. 

Since orthogonal polynomials of high orders are rapidly oscillating functions, a trun- 
cated expansion in conformal partial waves is, almost necessarily, oscillatory as well. Such 
a behaviour is clearly unphysical, but this does not constitute a real problem since physical 
observables are given by convolution integrals of distribution amplitudes with smooth coef- 
ficient functions. A classical example for this feature is the 77*-meson form factor, which is 
governed by the quantity 

where the coefficients are exactly the "reduced matrix elements" in the conformal expan- 
sion. The oscillating terms are averaged over and strongly suppressed. Stated otherwise: 
models of distribution amplitudes should generally be understood as distributions (in the 
mathematical sense). 

We give all relevant numerical input parameters for our model DAs in Table 01 at the 
scale fi = 1 GeV, which is appropriate for QCD sum-rule results, and, using the LO and NLO 
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0. 0.2 0.4 ,,0.6 



Figure 1: Left panel: 0| as a function of u for the central value of the hadronic parameters, for 
/i = 1 GeV. Red (solid) line: (p^.j^, green (long dashed): ^g.^, blue (short dashed): asymptotic 
DA. Right panel: same for ^g. 
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Figure 2: Left panel: 04;7r as a function of u for the central value of the hadronic parameters, 
for [I = 1 GeV. Red (solid) line: (p^.^, in conformal expansion, blue (dashed): 04;,^ using the 



renormalon model ^^t'^"'^ for the genuine twist-4 corrections. Right panel: same for 04;A'- 
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Figure 3: Same as Figure Elfor ip^. 
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scaling relations given in Section El and |31 at the scale fi = 2 GeV, in order to facilitate the 
comparison with future lattice determinations of these quantities. The mixing of i^T-meson 
parameters with operators of lower twist depending on is numerically small. 

The parameters related to twist-2 matrix elements have been determined using various 
methods; see the discussion in Section |21 Matrix elements of twist-3 and 4 operators for the 
TT meson were calculated a long time ago from QCD sum rules jEH 1121 HI] • In this paper, 
we perform a complete reanalysis of these sum rules and also include SU (3)-breaking effects 
relevant for the K meson. The corresponding sum rules and plots are given in the appendices. 
One important result is that we cannot confirm the sum rule for /s^ derived in Ref. and 
that our numerical value is considerably larger than that found in this paper. On the other 
hand, our central value for 6^ is similar to the one obtained in Ref. [51j, see also Ref. P7j . 

Finally, in Figure ^ we plot the twist-3 and -4 two-particle DAs for the vr meson, assuming 
massless quarks, and for the K meson, together with the corresponding asymptotic DAs. 
The figures show that quark-mass corrections significantly modify the end-point behaviour 
of 03, where they induce a logarithmic end-point divergency, even if the contributions of 
gluonic operators are neglected. This is not a problem because, as mentioned above, the 
DAs themselves need not be finite, it is only their convolution with perturbative scattering 
amplitudes that is meaningful. In Figures |21 and 01 we show the twist-4 two-particle DAs 04 
and also for the vr (left panels) and the K (right panels). The solid (red) curve in Figure|2l 
is obtained from Eq. ()4.32j) using the conformal expansion ()4.33|1 to NLO in the conformal 
spin, whereas the dashed (blue) curve includes the higher-spin contributions to the genuine 
twist-4 corrections as given by the renormalon model ()4.35|) . The mass corrections (p^^ 
vanish for the pion. It is clear that the higher-order contributions induced by ()4.35p modify 
both the end-point behaviour of 04;7r and the size of the DA away from the end-points. For 
the K, the absolute difference between both curves at, say, m = |, is very nearly the same 
as for 71, but the relative difference is much reduced because of large SU (3)-breaking effects 
induced by the mass-dependent contribution (p^x^ ■ Also note that the asymmetry of the 
curves induced by the non-vanishing value of af^ is not very pronounced, which is due to the 
smallness of that parameter as compared to = 1 and af". In Figure El we plot -04, with 
the same meaning of the curves as in Figure |21 Also here it is obvious that the renormalon 
model modifies the end-point behaviour of the DA, in particular for ip4-K, where it changes 
the sign of the logarithmic divergence at m = 0. 

6 Summary and Conclusions 

In this paper we have studied the twist-3 and -4 two- and three-particle distribution am- 
plitudes of ii'-mesons in QCD and expressed them in a model-independent way by a min- 
imal number of non-perturbative parameters. The work presented here is an extension of 
Refs. |17| 122 1^ and completes the analysis of SU (3)-breaking corrections by also including 
G-parity-breaking corrections in — rriq. Our approach consists of two components. One 
is the use of the QCD equations of motion, which allow dynamically dependent DAs to be 
expressed in terms of independent ones. The other ingredient is conformal expansion, which 
makes it possible to separate transverse and longitudinal variables in the wave functions, 
the former ones being governed by renormalization-group equations, the latter ones being 
described in terms of irreducible representations of the corresponding symmetry group. We 
have derived expressions for all twist-3 and -4 two- and three-particle distribution ampli- 
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tudes to next-to-leading order in the conformal expansion, including both chiral corrections 
0{ms + rrig) and G-parity-breaking corrections 0{ms — mg); the corresponding formulas are 
given in Sees. El and IH We have also generalized the renormalon model of Ref. |21^ to describe 
S't/(3)-breaking contributions to high-order conformal partial waves. 

We have done a complete reanalysis of the numerical values of the relevant higher-twist 
hadronic parameters from QCD sum rules. Our sum rules can be compared, in the chiral 
limit, with existing calculations for the vr |5H We confirm the sum rule for the twist-4 
matrix element 6^ quoted in Ref. [51], but obtain different results for the twist-3 matrix 
elements given in Ref. which lead to a 50% increase in the numerical value of the 
coupling /stt. Whenever possible, we have aimed at determining these matrix elements from 
more than one sum rule; we find mutually consistent results, which provides a consistency 
check of the approach. We have also studied the scale-dependence of all parameters to 
leading-logarithmic, or, if possible, next-to-leading-logarithmic accuracy, taking into account 
the mixing with operators depending on the strange-quark mass m^. Our final numerical 
results, at the scales 1 and 2 GeV, are collected in Table 01 

We hope that our results will contribute to a better understanding of S'[/(3)-breaking 
effects in hard exclusive processes and in particular in the decays of B and Bg mesons into 
final states containing K mesons. 
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Appendices 

A Non-Local Operator Identities 

For completeness, we quote the following non-local operator identities from Ref. |5Uj : 




d 




+ {mq - ms)q{x)ij5s{-x), 



(A.l) 





+ {rris + mq)q{x)i-f5s{-x) 



(A.2) 



df,q{x)ap^-f5s{-x) 



.8 



dx, 



V 



?(a;)75s(-x) + / dvvq{x)xpgGpy{vx)-ir,s{-x) 





+ {ms - mq)q{x)'yy-f5s{-x), 



(A.3) 
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q{x)a^u75s{-x) = -iduqix)'y5s{-x) + / dv q{x)xpgGpy{vx)'-f^s{-x) 



i / dvvq{x)xpgGpf,{vx)ap,y-f5s{-x) 



- {rus + mq)q{x)-iy'^^s{-x). 



(A.4) 



Here dp is the total derivative defined as 



dp {q{x)Ts{-x)} 



_d_ 



{q{x + y)[x + y, -X + y]Ts{-x + y)} 



By taking matrix elements of the above relations between the vacuum and the meson state, 
one obtains exact integral representations for those DAs that are not dynamically indepen- 
dent. 



B Sum Rules for Twist-2 Matrix Elements 



In this appendix we list and evaluate the QCD sum rules for twist-2 matrix elements of the K. 
The sum rule for fx, including 5'f/(3)-breaking corrections, was calculated in Refs. |^ ITU], 
that for af- in Refs. jHl CD] , and that for af- in Ref . [7] , apart from the perturbative terms 
in ml and the radiative corrections to the quark condensate, which are new. The sum rules 
read: 
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(gg) = (-0.24 ± 0.01)^ GeV^ 
{qagGq) = {qq} 


(ss) = (1 - 
(sagGs) = (1 - 


S3) {qq) 
d^){qagGq) 


(^^G^^ = (0.012 ± 0.006) GeV^ 






m2 = (0.8±0.1)GeV2, 5-^ 


= 0.2 ±0.2, §5 


= 0.2 ±0.2 


m,(2 GeV) - (100 ± 20) MeV ^ 


m,(lGeV) 


= (137 ±27) MeV 


as{mz) = 0.1187 ±0.002 


< — > as{l GeV) = 


- 

- U.JO_Q Q5 



Table A: Input parameters for sum rules at the renormalization scale = 1 GeV. The value of iris 
is obtained from unquenched lattice calculations with Nj = 2 flavours as summarized in Ref. |45] . 
which agrees with the results from QCD sum-rule calculations [47j. as{mz) is the PDG average 
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We evaluate the sum rules using the input given in Table El The results for fx and af^ are 
shown in Figure IU fx depends rather sensitively on the choice of sq. In order to reproduce 
the experimental result fx = 160 MeV, one has to choose sq = 1.1 GeV^. This is the value 
we will use also for all other sum rules for K matrix elements. For af^, we then find 



IGeV) = 0.30 ±0.15, 



(B.4) 



which is slightly larger than the result obtained in Ref. [7j and agrees with that obtained in 
Ref. jH]. For af- , we obtain the same result as Refs. [TIH ITT]: 



af (IGeV) = 0.06 ±0.03. 



(B.5) 



C Sum Rules for Twist-3 Matrix Elements 

In this appendix we estimate the parameters of the twist-3 distribution amplitudes f^K, ^3K 
and uj^K from QCD sum rules. Our approach is similar to that of Ref. |52j , where /s^r and uj^t^ 
have been determined, and based on the calculation of the correlation function of a non-local 
light-ray operator, which enters the definition of the three-particle distribution amplitude 
()3.3|1 . with the corresponding local operator: 

Un = t [ d%e-'Py{0\Tq{z)ta^,^,gG^,ivz)siO)siy)ta,,^59G.M'liyM 
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M 2 M 2 

Figure 4: Left panel: fx as function of the Borel parameter for sq = 1.1 GeV^. Solid line: 
central values of input parameters, dashed lines: variation of fx within the allowed range of input 
parameters. Figure taken from Ref. ^lOj. Right panel: same for af^ at the scale /i = 1 GeV. The 
results for af- are new. 



= {pz f J Vae-'P'^''^-^'"'-'^ TiDia) . (C.l) 
We also study the correlation function of that operator with the pseudoscalar current:^ 
IInd = ^ [d^ye-'Py{0\Tq{z)a,,^,gG^,{vz)s{Q)s{y)^,q{y)\Q) 



= {pzf j i)« e-'P"(°^+''"3) T^Noia) ; (C.2) 

for brevity, we do not show the Wilson lines in the non-local operators. Our calculation 
goes beyond that done in Ref. by including S'?7(3)-breaking corrections, and by also 
studying sum rules based on the non-diagonal correlation function, which allows a non-trivial 
consistency check of the results. 

Somewhat imprecisely, we will refer to Ud and Und SiS "diagonal" and "non-diagonal" 
correlation functions, respectively. The hadronic representation of the non-diagonal correla- 
tion function Und only contains pseudoscalar = 0~ contributions, whereas the diagonal 
correlation function Ho also contains contributions of states with higher spin, = 2~ and 
J'^ = 1+. This is not a disadvantage, since such states all have considerably higher masses 
than the K meson, and can effectively be thought of as parts of the continuum contribu- 
tion. For reasons that will become clear below, we have also calculated a correlation function 
similar to ()C.1|) . but with operators of opposite parity: 



UD = t J d^ye-'Py{0\Tq{z)a^,gG^,{vz)s{0)s{y)a,,gG,M(l{yM 

= {pzf j Pa e-'P"("^+""^) T^D{a) . (C.3) 



®Note that the currents in IIjvd contam no factors i, in contrast to Ud. This is so as to obtain a positive 
spectral density. 
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For the diagonal correlation function we find, using factorization approximation for the four- 
quark condensates and dropping terms that vanish after Borel transformation: 
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To this accuracy, the expressions for ttd and ttd are almost the same, the only difference 
being that in vfo the last term in ()C.4jl . the contribution of {qq){ss), comes with a minus 
sign. In the chiral limit, we can compare the above result with that obtained in Ref. [S2]: we 
find agreement for the perturbative contribution, but a different answer for the contribution 
of the four-quark condensates. The leading-order contribution 0{as) of the gluon condensate 
as well as that of the dimension-6 triple-gluon condensate (g^fG^) both vanish. We also have 
calculated the contribution of the gluon condensate in the local limit, Q-^p^i'^^+vas) _^ ^^^^ 
find 
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which differs from the result obtained in Ref. In particular, we do not reproduce the 

logarithmic term quoted in |5"2j . 

For the non-diagonal correlation function we find 
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The sum rules for the couphngs /sx, As^- and u^k are derived by expanding the correlation 
functions in powers of {pz): 

Un = {pzf {ng) + i{pz) [ng) + {2v - i)ng)] + o{{pzf)] , 
HiVD = (p;^)' {ngg + ^(pz) [n!Jg + {2v - i)n!;;g] + o((pz)^)} . (c.7) 

Comparing these expressions with the corresponding expansion of the K contribution to the 
correlation functions expressed in terms of the DA ()3.5|) . one obtains 
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from the diagonal and non-diagonal correlation functions, respectively. Here and below 
B[. . •](M^) stands for the Borel transformation with respect top^; is the Borel parameter. 
From II/), we obtain the following sum rule for f^x'- 
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where the subscript D indicates that this sum rule is derived from the correlation function 
Hd. The last term on the right-hand side comes from the factorisation of the four-quark 
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condensate {qcr^ut^q){qafj^iyt'^q). In Ref. the authors have argued that this term, which 
induces a large power correction in their sum rule for /a^, is unreliable because of a potential 
breakdown of the factorisation approximation for that particular condensate; they suggested 
to determine /s^ from a sum rule derived from the sum of the correlation functions Ud + TIo 
instead, where these large contributions cancel. Indeed, the Dirac structures af^u and ia^y-^^^ 
are not independent, but related by icr^^'j^ = —^e^upaO'pa, which induces the relation 

^D = tJ d%e-'^y{0\Tm^a^zl69G,MsiO)siy)^cT^zl,9G,,iy)qiy)\0) 

- I j Ae-*P^(0|Tg(0)za^,75^7G',.(0)s(0)s(i/)za,,75(7G'^.(2/)g(l/)|0) . 

IId receives no contributions from 0^ states because their contributions to the two correlation 
functions on the right-hand side are equal and cancel in the difference; the same applies to 
1+ states, so that the lowest resonance contributing to lir) is 1~. These states can safely be 
included in the continuum so that it is possible to extract f^x from the sum of correlation 
functions TId + ^d- On the other hand, our sum rule ()C.10|) . derived from only, with the 
correct coefficients for gluon and four-quark condensates, is actually not very sensitive to the 
term in {qq){ss), but dominated by the gluon condensate. As there is no strong theoretical 
argument in favour or disfavour of either diagonal sum rule, the one based on II^) and the 
one based on Ud + Tljy, we decide to use both. We also determine f^K from a third sum 
rule based on the non-diagonal correlation function IliVD; the difference between these three 
results will be interpreted as theoretical uncertainty. 

Explicitly, we obtain, in addition to ()C.10|) . the following sum rules for f^x, with the 
index indicating the underlying correlation function: 

89 /«« 1 



2 a., , > , /'S , r da 
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+ -^^i.^s{s(TgGs) + mg{qagGq)) 
The sum rules for f^^^ are obtained by taking the chiral limit of the above expressions. 
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As for uJsK, we have not calculated the gluon-condensate contribution to the diagonal sum 
rule, which is expected to be dominant, so we cannot use the diagonal sum rule and only 
consider the non-diagonal one: 



nis + mq oOtt'^ Jo 27 tt 
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In evaluating this sum rule, we replace fj,K by the expression obtained from ()(I12jl . 

As for \-iK-, the gluon-condensate contribution is suppressed by a factor — by 
virtue of G-parity and can safely be neglected in the diagonal sum rule. We did calculate 
this contribution for the non-diagonal sum rule, though, where indeed it gives only a small 
contribution. On the other hand, the (gg) (ss) contribution is also absent because of G-parity, 
so that the two diagonal sum rules for /l/^Aai^-l^ and /I^^Aai^l^^^j differ by a global factor 2. 
As the values of /l^-l^) and f^Klo+D ^^^^ differ by a factor of approximately 2, this theoretical 
uncertainty cancels to a large extent. The sum rules read: 
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Figure 5: Left panel: /s^r as a function of the Borel parameter, calculated from the non-diagonal 
sum rule (IC.12j) (red, solid line), the pure-parity diagonal sum rule (IC.lQj) (green, long dashes) 
and the mixed-parity diagonal sum rule (IC.llj) (blue, short dashes); sq = 0.8 GeV^. Right panel: 
same for f^x', sq = 1.1 GeV^. 




Figure 6: Left panel: lo^j, as a function of the Borel parameter from the non-diagonal sum rule 
(f03|l : So = 0.8 GeV^ Right panel: same for ujsk', sq = 1.1 GeV^. The results from diagonal 
sum rules are not shown because the gluon-condensate contribution is unknown. 
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Figure 7: X^k as a function of the Borel parameter, calculated from the non-diagonal sum rule 
()C.16^1 (red, solid line), the pure-parity diagonal sum rule (|C.14p (green, long dashes) and the 
mixed-parity diagonal sum rule ()C.15|) (blue, short dashes); sq = 1.1 GeV^; As^ = by virtue of 
G-parity. 
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Again, when evaluating these sum rules, we replace f^K by the corresponding expressions 
obtained from ^J^, KHA^ . and (im2|l . 

The numerical results from all these sum rules are shown in Figures El to [3 As for /s^r 
and fsK, all three sum rules yield very similar results, which is a strong indication for the 
consistency of the approach. The diagonal sum rules are very stable in M^, the non-diagonal 
ones less so. Taking into account the uncertainties of the input parameters as given in Table 1X1 
and the difference in the results from the different sum rules, we obtain the estimates 

/3^(1 GeV) = (0.0045 ± 0.0015) GeV^ f^Kil GeV) = (0.0045 ± 0.0015) GeVl (C.17) 

The effect of SU (3) breaking is very small, 

W/3. = 0.98 ±0.03, (C.18) 

as all sum rules are dominated by the contribution of the gluon condensate. Note that our 
value for /s^ is about 50% larger than the one obtained in Ref. [H2], which is due to, as we 
believe, the incorrect results for the contributions of the gluon and four-quark condensate 
contributions obtained in this paper. As for 003, as explained above, we only evaluate the 
non-diagonal sum rule. We find that the sum rules are less stable in M^, as with /a, and 
that now the effect of SU{3) breaking is more prominent. Our final estimate is 

c<;3^(lGeV) = -1.5±0.7, o;3i^(l GeV) = -1.2 ± 0.7 , (C.19) 

where the error reflects in particular the uncertainty of the value of the gluon condensate. 
Our result is to be compared with that of Ref. cua^ ~ —3. Finally, X^k can be determined 
from three sum rules, as the gluon-condensate contribution is suppressed by a factor m^. All 
three sum rules yield perfectly consistent values, despite the fact that the two diagonal sum 
rules flG.14|) and ()G.15|) differ by an overall factor of 2, which, as expected, is largely cancelled 
by the different values of f^^lo and /I^Id+d- We obtain 

A3K(lGeV) = 1.6±0.4; (C.20) 

the error is smaller than for u^k because the gluon condensate is suppressed. This result is 
new. 

D Sum Rules for Twist-4 Matrix Elements 

The aim of this section is to estimate the decay constant Sj^ that determines the normalization 
of twist-4 distribution amplitudes. To this end we define the currents 

= qgG^.aiaS , Jl = qgG^^alalbS , (D.l) 
with quantum numbers = 1+ and 1~, respectively, and calculate the correlation functions 

K'^ = ^ j d'xe^^^ {^\TJ^y{x){J^y)\m = P.P. n^''^(p2) - nf '^(p^) , (D.2) 

taking into account contributions of operators with dimension up to eight. Note that the 
relative sign between fx and 5"]^ can be fixed from the non-diagonal correlation function of 
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and the axial vector current. This calculation was done in Ref. [21] and will not be repeated 
here; the result is that Sj^ is positive. 

Similar correlation functions have been considered in the past, mainly in connection with 
searches for exotic quark-ant iquark-gluon mesons We obtain 
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In both cases the upper sign refers to the axial and the lower sign to the vector correlation 
function, respectively; 11^ has been calculated, in the chiral limit, in Ref. jS^- The quark 
mass corrections and the expression for 11^ are new. 
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Figure 8: Left panel: 5^ as a function of the Borel parameter from the sum rule (|D.5p : sq = 
0.8 GeV^ Right panel: the same for 5l; Sq = 1.1 GeV^ 

In this work we follow the procedure proposed in Ref. jM] and write the sum rule directly 
for the correlation function 11^: 

The results for and S]^ are shown in Figure |H1 We find 

6l = (0. 18 ± 0.06) GeV^, 5| = (0.20 ± 0.06) GeV^ (D.6) 

and 

4/5^ = 1.10 ±0.05, (D.7) 
which can be compared to the estimate (/E'<5|-)/(/7r5^) = l-07lo;}3 obtained in Ref. |23j. 
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